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Interpreting a morphogen gradient into a single stripe of gene-expression is a fundamental unit of 
patterning in early embryogenesis. From both experimental data and computational studies the 
feed-forward motifs stand out as minimal networks capable of this patterning function. Positive feedback 
within gene networks has been hypothesised to enhance the sharpness and precision of gene-expression 
borders, however a systematic analysis has not yet been reported. Here we set out to assess this hypothesis, 
and find an unexpected result. The addition of positive-feedback can have different effects on two different 
designs of feed-forward motif- it increases the parametric robustness of one design, while being neutral or 
detrimental to the other. These results shed light on the abundance of the former motif and especially of 
mutual-inhibition positive feedback in developmental networks. 



Modularity is a central principle to both systems and synthetic biology 1,2 . It implies the reduction of large 
gene networks to a collection of small, separable subsystems. Milo et al. 3 initiated such a collection by 
identifying small groups of interacting genes or network motifs that appear statistically overrepresented 
in complex biological networks. This statistical motif-extraction usually involves the structure - but not the 
function - of the motif. A proper modularity requires that these motifs always perform a similar function in 
different biological contexts 4 6 . For such an assessment, systems and synthetic biology increasingly aim to explore 
the biological and dynamical functions of these minimal regulatory motifs. Existing studies mostly addressed the 
characteristics and requirements of genetic oscillators, biological switches and low-pass filters 7 9 . As the list of 
well-studied motifs is growing, several questions emerge: What happens when simple modules are combined 
together? Does a given motif perform the same function, whether in isolation or embedded in a large network? 
(Fig. 1). 

Here we address these questions by studying the motifs capable of performing a specific biological task. The 
task chosen is converting a decreasing morphogen gradient into a single stripe of gene expression. The single 
stripe is a key spatial feature in the early stages of the embryogenesis. It is considered a mechanistic unit in the 
formation and maintenance of boundaries between future tissues 10 . This scenario of an input gradient determin- 
ing distinct cellular fates is referred to as positional information 11 . As a definition, a signalling molecule that 
generates at least two distinct cell types, or thresholds, at different concentrations is called a morphogen 12 . 
Examples of such threshold-controlling morphogens include Bicoid and Dorsal in Drosophila u,M , activin in 
Xenopus 15 and Sonic Hedgehog in the vertebrate neural tube 16 . 

Initially, it was hypothesised that these thresholds are determined by different responsiveness (i.e. different 
number and affinity of binding sites) of the genes to the morphogen. However, recent evidence suggests that it is 
the regulatory logic of the developmental network itself that determines the expression patterns 14,17 " 20 . In other 
words, the response to morphogen gradients is an emergent property of the regulatory circuit rather than a 
property of specific genes. Together with the identification of spatiotemporal motifs from developmental gene 
networks 21 , these observations suggest that general circuit-design principles are responsible for pattern formation 
at tissue level. Studying simple but essential expression patterns and the genetic circuits that produce them 
constitutes an approach to uncovering such principles. 

From the existing studies of gene networks capable of a single-stripe pattern under a morphogen gradient, the 
3-gene Incoherent Feed-Forward Motifs (IFFMs) appear as the minimal ones. These studies consisted of general 
computational explorations 22 24 or specific experimental implementations in synthetic and developmental bio- 
logy 2527 . More exhaustively, Cotterell & Sharpe 28 computationally analysed all 3-gene configurations and 
revealed that indeed the networks showing the single-stripe expression include the 3-gene IFFMs. 
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Figure 1 | Examples of functional motifs in biological networks. After identifying the biological and dynamical functions of individual network motifs, 
either genetic or biochemical 8 , a major question still remains: will a given motif have consistent effects when added to different networks? 



Beyond the mere existence of a single- stripe pattern, the sharpness 
and precision of this simple gene-expression pattern are crucial for 
the correct cell-fate specification. An unmet challenge for devel- 
opmental biology is to elucidate the mechanisms that control these 
features and keep cell populations separated 29 . Among possible 
mechanisms, positive feedback within the gene network has been 
repeatedly hypothesised to be the primary mechanism 30 32 . Here, 
we set out to examine whether this hypothesis is correct, and thus 
assess the extent to which the positive feedback controls the forma- 
tion and features of the stripe pattern. Firstly, we analytically deter- 
mine the requirements and characteristics of stripe formation in 
IFFMs, with emphasis on developmental gene networks. Secondly, 
our computational exploration shows that, in spite of their similar- 
ities, the motifs respond differently to the addition of positive feed- 
backs. Finally, we explain this unexpected behaviour by coupling 
network architecture and dynamics. 

Results 

The model. We employ a model in which a one-dimensional tissue 
receives a morphogen gradient that is transformed by the cells' gene- 
tic network into a spatial pattern, here a single stripe expression in 
one of the genes (Fig. 2). Represented as a square, gene A is the only of 
the three genes to receive the input signal. Gene B is considered 
throughout this study as the stripe-gene, while gene C is the 
intermediate gene. 

Genes integrate multiple regulations from other genes following 
the connectionist model 33 . This model is biologically-verified and 
widely- employed, especially in studying developmental gene net- 
works 13,28,34 ~ 36 and signalling networks 8 . The regulation from gene 
X to gene Y is characterised by a sign (activating or inhibiting) and 
a weight (co X y) (see inset from Fig. 2). In the connectionist model, the 
effective regulating input to a gene is the sum of the weighted con- 
tributions. Finally, the rate of protein production from this gene is 
proportional to a saturating function of this input. We refer to this 
approach of summing the weighted contributions and applying a 
saturating filter- function as the sum & filter regulation model. 

In endogenous tissues, proteins can diffuse between cells. How- 
ever, Jaeger et al. 17 found that diffusion is not required for qualita- 
tively generating the expression pattern of the gap genes in 
Drosophila. Cotterell & Sharpe 28 reached the same conclusion by 
means of an exhaustive computational study of all stripe-forming 
3-gene networks. Therefore, the current study considers that 



proteins are non-diffusible, making our framework a positional- 
information scenario. 

Two motifs constitute the minimal band-generators. There are 
four possible IFFM configurations, denoted as 1 1 - 1/ (Fig. 3). 
Below each of the four minimal topologies, graphs illustrate how 
each motif can, or cannot, make a stripe. The upper graphs include 
the individual contributions from genes A and C to gene B, under 
increasing morphogen input, M. The sum of the two contributions 
constitutes the effective input into gene B. This effective input is 
transformed into the output (lower graphs) by the sigmoidal 
regulation function and time integration (see Methods section). 

These graphs reveal two main aspects of the signal integration in 
IFFMs. Firstly, symmetry is observed in the signal integration for the 
stripe-forming T and I 3 motifs: the two contributions to B are sym- 
metric with respect to the morphogen axis. The symmetry stems 
from the topology of the two paths from morphogen to stripe gene: 
direct positive & indirect negative in I u and indirect positive and 
direct negative in I 3 (see Supporting Information, Fig. SI). 

Secondly, the signal integration for the I 2 and I 4 motifs sheds light 
on the patterning capabilities of these motifs. I 2 has only negative 
interactions. According to the connectionnist model, these negative 
inputs are transformed into low activity from the promoter, as the 
promoters are considered to be non-constitutive. In other words, the 
I 2 motif is able to show at most a barely-detectable stripe. 

A few comments on the behavior of I 2 are in order. The current 
study builds upon the results and accordingly borrows the connec- 
tionist framework of Cotterell & Sharpe 28 who identified the stripe- 
forming 3-gene networks and classified them into categories 
characterized by a common dynamical mechanism. In spite of the 
generalities of the model, a few assumptions are worth noting for 
their effect upon the conclusions. In the absence of external activa- 
tors, self-activation can compensate for the lack of constitutive genes. 
From this perspective, a non-constitutive 3 -link I 2 network can 
not produce a stripe in the absence of self-activations. In other 
words, a constitutive 3-link I 2 or equivalently a non-constitutive 5- 
link I 2 (Fig. 3) can show a pronounced single stripe of gene 
expression. 

Finally, among the four IFFMs, the I 4 motif stands out as the only 
motif that cannot show a stripe under the current non-constitutive 
sum & filter approach. We qualitatively show this feature in Figure 3, 
while Supporting Information includes the detailed analytical proof. 
From these results together with the observation of low output from 



SCIENTIFIC REPORTS | 4:5003 | DOI: 1 0. 1 038/srep05003 



2 



Exponential 
Morphogen 
Gradient 



Tissue of 
N Isogenic 
Cells 




Regulation Function 



Cell Response 

through 

Gene 

Expression 



• B 

' C ( 









Cell index 



Interaction matrix: 





M 


A 


B 


c 


M 


0 


0 


0 


0 


A 


1 


0 


0 


0 


B 


0 




0 


UCB 


C 


0 


<+>AC 


0 


0 



Synthesis rate of the protein 
is proportional to: 



/(*) = 



1 

1 + exp(a — b x) 



Additive input into e.g. B: 

input, x x = ojabA + u)cbC 



Figure 2 | The framework of the current study. An example of single-stripe formation in the positional-information scenario: the concentration of a 
morphogen Mis monotonically increasing along a one-dimensional "tissue" of N isogenic cells. The underlying genetic network transforms the 
increasing input into a single stripe in the expression of one of the genes. Each transcription factor X acting at the promoter region of gene Yhas an effect 
on gene expression characterised by the interaction strength co X y (Interaction matrix). The activating interactions have oj X y > 0 and inhibiting 
interactions, oj X y < 0. For each gene, the contributions from multiple transcription factors affecting it are subsequently summed. The resultant 
transcription rate from the promoter is proportional to the sigmoidal-filtering of this total contribution. The parameters a and b of the sigmoidal function 
fx) control the steepness and location of the threshold value of the regulation function. The morphogen Mis considered to be received only by gene A (see 
also Methods section). 



I 2 network, we conclude that only and I 3 are core topologies for 
single-stripe formation. 

Quantitatively, we also obtained analytical approximations of the 
factors controlling the location and width of the single stripe in 
IFFMs (Supporting Information, Fig. S2). Firstly, the existence of 
the peak is conditioned by a balance between the strengths of the 
activating and the inhibiting path into gene B. The required balance 
of weights is modulated by the protein degradation rate and the 
steepness of the regulation function (Eq. (7) in the electronic supple- 
mentary material). Secondly, the location of the peak is uniquely 
controlled by the weight of the intermediate interaction, oj ac . This 
weight is drastically limited to low values. This ensures the position- 
ing of the stripe within the one-dimensional tissue, as high values of 
co AC locate the peak outside the posterior end of the 'tissue'. 
Additionally, the positive contribution to B has an upper limit as 
well, as it determines the width of the stripe. It is co AB in l l and co CB in 
I 3 . In both motifs, the negative interaction (co C b in Ii and co AB in I 3 ) is 
not intrinsically restricted, being limited only with respect to the 
positive interactions. This relative restriction stems from the sum 
& filter regulation: promoters respond actively only in the presence 
of activators, whose effect is diminished or entirely cancelled by 
inhibitors. Numerical simulations reflect that indeed the positive 
interactions are limited to low parameter values while the negative 
ones have no upper limit (Supporting Information, Fig. S3). 

Thus, we have found the parametric requirements for the two 
minimal gene-networks to transform a monotonically distributed 
input signal into a band-like expression in one of its constituent 
genes. These motifs achieve the patterning task in a symmetric and 
similar manner. Subsequently, we inquired into the effect of positive 
feedback in these alternative designs of stripe formation. 

Constructing a topology-atlas by adding positive feedback. To 

assess the effect of positive feedback on the stripe formation, we 
studied a class of motifs that we refer to as RIFFMs: Reinforced 



Incoherent Feed-Forward Motifs. This collection (Fig. 3) consti- 
tutes the motifs based on IFFMs with all possible combinations of 
positive feedbacks added on and between the two downstream genes 
B and C. No feedback to gene A in these motifs would conserve the 
feed-forward nature of the signal. 

The representation of functional Ij- and I 3 -based RIFFMs (Fig. 3) 
recalls two 3-dimensional cubes whose nodes are the RIFFMs, and 
whose edges are the topological changes generating positive feed- 
backs. There are three types of edges associated to the addition of 
positive feedback. Two of them are direct positive-feedbacks: self- 
activation of gene B (dotted edges) and of gene C (solid edges). The 
third type is an indirect positive feedback formed by adding B-C 
interaction (dashed edges). In the case of 1^ this must be a double- 
negative, while in I 3 , it is a double-positive. 

To inquire which among the added interactions are beneficial, 
detrimental or neutral with respect to the interpretation of the mor- 
phogen gradient into a single stripe, we numerically determined the 
functional parameter space and the characteristics of the resultant 
stripes for all RIFFMs. Details on the definition of stripe are found in 
Methods. Figure 4 summarises the effects in terms of parameter- 
space extension brought by the additional interactions. This figure 
is a schematic, qualitative representation of the full quantitative, 
numerical calculations of all motifs included in Figures S4, S5 of 
Supporting Information. The graphs in Figure 4 show the compar- 
isons performed between neighbouring RIFFMs (i.e. networks con- 
nected by one edge in the cube-representation from Fig. 3). For each 
branch (Ij and I 3 ) we have illustrated the methodology of this com- 
parison by chosing one of the three possible new feedbacks, and 
indicating (qualitatively) how the size and position of the functional 
parameter space (red region) changes as the new regulatory link is 
gradually added to the network (lower part of 4). In both cases we 
have chosen the indirect feedback case - ie. the feedback between B 
and C (the mutual repression for 1 1} and the mutual activation for I 3 ). 
These examples show that the functional region in (oj ab , (o AC , Mcb) 
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Figure 3 | The formation of the single-stripe pattern by 3-gene networks based on the four Incoherent Feed-Forward Motifs. The upper part of the 
figure shows the stripe-forming IFFMs reinforced with positive feedbacks, having from 4 to 6 interactions. Only interactions that generate positive 
feedbacks have been included. Within the i! and I 3 branches of stripe-forming networks, thick solid-line arrows imply adding C self- activation, the 
dashed-line arrows, B self-activation, and the dotted-line arrows, mutual-interaction (i.e. mutual-inhibition or mutual-activation). The lower rectangle 
in the figure includes the four IFFMs and the underlying process of stripe formation. It shows the generic stripes of the Ii and I 3 networks, the lack of stripe 
for I 4 and the barely-detectable one for I 2 . Within the rectangle, the upper graphs single out the individual contributions into the stripe-gene promoter for 
an increasing morphogen, M: from gene A in red, and from gene C in blue. For simplicity, we chose to use the concentration of A for the horizontal axis, 
instead the morphogen concentration. The lower graphs show the resultant expression from gene B. The Gray areas represent non-biological negative 
values for the concentration A and conceptually illustrate that I 4 would require a negative concentration of A to show a stripe. The range of (non- 
biological) negative concentration for protein A is also included for illustrative purposes. 



expands when adding the B-C inhibition to I,, and remains comple- 
tely unchanged when adding B-C activation to I 3 . In other words, the 
latter behavior implies that any functional set of (co AB , m AC , (o C b) will 
maintain its functionality when the B-C activation is added. The 
examples in the lower part of Figure 4 compare only two 3-link 
networks to one of their 4-link neighbours. The full analysis schem- 



atically summarised in the upper graphs of Figure 4 consists in the 
comparison of the functional parameter space of all relevant motifs 
from one level of the atlas (i.e. w-link networks, with n being 3, 4 or 5) 
to the level above (i.e. (n + l)-link networks). 

We found that within each branch all these comparisons qualita- 
tively group into three classes that directly relate to the three types of 
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Figure 4 | Changes in the functional parameter space produced by adding positive feedbacks to the IFFMs. The positive feedbacks added to the IFFMS 
are: two self-activations (of gene B and C), and the B-C interaction. We show thus three graphs associated to these three added positive feedbacks. The 
graphs illustrate the change in the extent of the functional parameter space as the positive feedback increases. More precisely, they include the fraction of 
the functional parameter space of the RIFFMs, each scaled to the extension of the functional parameter space of a lower level RIFFM. This comparison 
methodology is illustrated in the lower panels for the mutual-inhibition and the mutual-activation of the 4-link networks. In these two examples, the 
parameter space to be compared is (co^s, tojo cocb) with (absolute) values co X y e [0, 10] (Gray cube). The functional region of parameter space is 
represented by the red shape. The behaviour illustrated in these three graphs is generic for the addition of each of the three types of positive feedback, 
independent of the existing positive feedbacks in the network. We wish to point out the parametric boost brought by the mutual-inhibition in the Ij 
branch, and the complete neutrality introduced by the mutual-activation and stripe-gene self-activation in the I 3 branch. 



added interactions: B self- activation and C self-activation, and the 
B-C interaction (Supporting Information, Fig. S4, S5). The qual- 
itative classification of behaviours is illustrated in the three upper 
graphs in Figure 4. This qualitative grouping of comparisons into 
three classes implies that the addition of an interaction generating 
positive feedback has a similar effect when introduced either to a 3- 
link, to a 4-link or to a 5-link network, in the same branch. In other 
words, for a given motif, the effect of adding positive-feedback is 
independent on the existing interactions. We thus say the positive 
feedback is modular within a motif. In contrast to this similarity 
within a motif, Figure 4 reflects that positive feedback has fun- 
damentally distinct effects in the two branches, Ij and I 3 . We thus 
say that the same positive feedback is non-modular between differ- 
ent motifs. In the following we discuss in more details these non- 
modular differences. 

Parametric-boost of the I t design upon addition of positive feed- 
back. In the Ii branch, as the strength of each additional feedback 
interaction is gradually increased, the functional parameter space 
firts grows in size until it reaches a maximum and then declines. 
Regarding the pattern itself, within the entire branch of networks, 
introducing and increasing B self- activation widens and sharpens the 
single stripe pattern. Introducing and increasing C self-activation has 
an effect primarily on the right border of the stripe, displacing 
it posteriorly. In other words, strong values for self-activations 
shrink the functional parameter space. 



Most interestingly, the highest parametric boost was found to 
occur for the creation of the mutual-inhibition module, i.e. by the 
addition of the B-C inhibition. This boost is observed at any level of 
the topology atlas (i.e. it is equally observed whether or not the B and 
C self- activations exist). In other words, for any motif in the 1 1 branch 
with less than 6 interactions, the addition of the B-C inhibition 
enriches the functional space with parameter sets that were not func- 
tional in the absence of this interaction. While a small amount of 
functional parameters is lost with the addition of this interaction, 
significantly more parameter sets are added, producing thus a net 
gain of functional parameter space. Effectively, the functional region 
shifts in parameter space, as well as growing. 

The decrease in the fraction of parametric volume for the Ij branch 
beyond | co^c | ~ 4 in Figure 4 is a consequence of the limitation of the 
parameter space to the arbitrarily- chosen interval iO X y e [ — 10, 10] . 
For \o)bc\ ~ 4, the stripe formation requires positive-interaction 
values higher than the co Max =10 considered in the numerical study. 
Consequently, this decrease in the functional parameter space is not 
an intrinsic stripe-formation restriction (Supporting Information, 
Fig. S6). This observation emphasises the fact that indeed mutual- 
inhibition intrinsically expands the functional parameter space at 
any level (i.e. to a 3-link, 4-link or 5-link motif) it is introduced within 
the Ii branch. 

Parametric-invariance of the I 3 design upon addition of positive 
feedback. Fundamentally different from the I 1 branch, the I 3 branch 



SCIENTIFIC REPORTS | 4:5003 | DOI: 1 0. 1 038/srep05003 



5 



shows complete neutrality for all interactions except the self- 
activation of (intermediate) gene C, for which it shows fragility 
(Fig. 4). By neutrality we mean that parameters sets which are 
functional in the absence of B self- activation and/or B-C activation 
will remain functional when these interactions are added. Added 
separately or simultaneously, these two interactions widen the exist- 
ing stripe and sharpen its borders. Consequently, the functional 
parameter space of I 3 does not expand under the addition of 
positive feedback (see also Supporting Information, Fig. S7). 

As in the case of Ii branch, the addition of each interaction sepa- 
rately is modular within the branch. For example, the addition of C 
self-activation shows the same detrimental effect within the entire 
branch (i.e. either introduced at level 3 or at level 5), following the 
behaviour illustrated in the upper right corner of Figure 4. 

Network topology drives the dynamical response to node- 
reinforcement. There are clear similarities between the Ij and the 
I 3 : they interpret the morphogen gradient into a single-stripe pattern 
under similar parameter conditions, display similar topological 
complexity, and operate in a symmetrical manner. Given these simi- 
larities, one might have expected modularity between these two 
branches (in addition to the modularity seen within the branches). 
Consequently, we sought a mechanistic explanation of these intri- 
guing results, mainly of the parametric boost in the I] branch and the 
neutrality in the I 3 branch. 

We found the explanation to be twofold. Firstly, the parametric 
boost brought by B-C inhibition in T branch can be traced back to 
the sum & filter modelling approach. We have established above that 
a strong restriction exists for oj ac in functional Ij IFFM. High values 
of (Oac push the stripe's center outside the limits of the 'tissue', 
resulting into the loss of stripe. When mutual-inhibition exists, C 
receives an activation and an inhibition that can compensate each 
other and achieve the functional input value to gene C (i.e. to AC A + 
oj bc B). In other words, iO AC can now acquire a wider set of values, as 
it will be compensated or reduced by B-C inhibition. This wider set of 
values is responsible for the parametric boost produced by the addi- 
tion of B-C inhibition. 

However, this compensating- interactions reasoning does not 
account for the neutrality present in the I 3 branch. On the contrary, 
fragility would be expected instead. In consequence and as the second 
part of the mechanistic explanation, we analysed the observed beha- 
viour from the dynamic systems perspective. In the following, we 
briefly define several useful concepts from dynamic-systems theory, 
such as attractors, phase-portraits and nullclines. 

The expectation that single-stripe formation should benefit from 
the addition of positive feedbacks stems from the conferred multi- 
stability, i.e. co-existence of multiple stable steady states 30,31 . A steady 
state is a specific value of the system's variables (here the concentra- 
tions A, B, C) that is indefinitely maintained by the system, in the 
absence of perturbations. Such a state is said to be stable if, under 
small perturbations, the system returns to this steady state; and it said 
to be unstable, otherwise. The stable steady states, or in short attrac- 
tors, define distinct cell fates, a view proposed more than 60 years ago 
that has received recent confirmation 37 . In a multistable system, the 
switching from one attractor to another can act as an effective mech- 
anism for sharp-boundary generation in the early developmental 
stages 30 . The attractor- switching process in itself generates two 
regions with distinct identities, conferred by distinct attractors, con- 
stituting the ideal mechanism for the separation of body segments 
and tissues. 

In the current study, the lack of positive feedback in the four 
IFFMs implies the existence of a sole attractor (Eq. (1) in Support- 
ing Information). Therefore, the stripe shape is the result of the 
attractor increasing and subsequently decreasing its _B-value in res- 
ponse to the input signal. In other words, IFFMs' stripe has borders 
due exclusively to the movement of this sole attractor in phase space. 



The borders of this stripe are generically smooth and shallow. By 
containing positive feedbacks, RIFFMs can show multistability. In 
this case, in addition to moving the steady states, the change in 
morphogen concentration along the tissue also results in the destruc- 
tion or creation of new steady states. By this change, RIFFMs show a 
qualitatively new type of stripe for which one or both borders are 
formed by switching between two co-existing attractors. As a repres- 
entative example, the study by Manu et al. 36 concerning the stripe 
formation by the gap genes in Drosophila is the first to make the 
distinction between gene-expression borders formed by attractor- 
movement or by attractor-switching. 

To exemplify these concepts, Figure 5 shows the generic stripe 
produced by adding mutual-interaction to I x and I 3 . A graphical tool 
for determining the steady states is drawing the system's nullclines in 
the state space, here the space of the protein concentrations B and C. 
The nullclines are defined as the curves whose points are charac- 
terised by one variable not changing in time, i.e. the B-nullcline 
contains all those points for which the system will maintain B con- 
stant. One can have an intuitive understanding of the dynamics of the 
system from looking at the 2-dimensional (£, C) subspace instead of 
the full 3-dimensional (A, B, C) space. In such a simplified visualisa- 
tion (upper and lower graphs in Fig. 5), A is considered to be at steady 
state, rapidly following the morphogen gradient. This simplification 
is allowed by the feed-forward nature of the RIFFMs. This 2D sim- 
plification only serves for a qualitative understanding of the typical 
behaviour of a RIFFM, as the timecourse of a specific parameter set 
requires the study of the 3D system. Employing this 2D space, we 
studied the change in shape and position of the nullclines produced 
by the change in the morphogen concentration, and the conse- 
quences on the associated steady states and thus stripe formation 
(see also Supporting Information, Movies M1-M4). 

We found the shape of the nullclines to be determined exclusively 
by the interaction strengths, a) XY , while the position of the nullcline in 
the (C, B) space is essentially controlled by the morphogen gradient 
(electronic supplementary material). More precisely, the change in 
morphogen concentration leads to a simple horizontal or vertical 
shift of the nullclines, as depicted by the red arrows in the anterior 
(C, B) graphs in Figure 5. The shift of nullclines creates and annihi- 
lates steady states. The simple nature of this shift provides a mech- 
anistic perspective on the creation, destruction and selection of 
attractors, and the use of these events by the stripe-formation 
process. 

The examples in Figure 5 constitute the generic response of the two 
motifs, Ij and I 3 , to the addition of mutual- interaction positive feed- 
back. There, a non-functional parameter set for T IFFM becomes 
functional with the addition of B-C inhibition, leading to a sharp 
band (middle panels). Instead, the addition of B-C activation to the I 3 
IFFM sharpens a stripe already present in the IFFM. 

Finally, we turn to the explanation of the neutrality in the I 3 
branch. As exemplified in Figure 5, the addition of B-C activation 
to the I 3 IFFM changes the shape of the blue C-nullcline from a 
simple vertical curve to a sigmoidal. This transformation allows for 
multiple steady states, with the attractors generally located on the / 
diagonal in the (C, B) space. The stronger the interaction, the steeper 
the C-nullcline. Nevertheless, this steepness simply translates into a 
rapid attractor- movement for the left border, and a sharp attractor- 
switching for the right border of the stripe. In other words, this 
configuration of attractors allows an unlimited increase in the B-C 
interaction. Similarly, Supporting Information Figure S8 shows how 
the nullclines' configuration allows for full neutrality of the single- 
stripe pattern to the addition of B self-activation. 

In the Ii branch as well, the addition of B-C inhibition changes the 
shape of the blue C-nullcline from a simple vertical curve to a sig- 
moidal. However, the attractors are now located on the \ diagonal in 
the (C, B) space, instead of the / diagonal. This nullclines' configura- 
tion favours the attractor-switching for both borders. Small values of 
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Figure 5 | Parametric boost versus neutrality: a mechanistic view. The figure illustrates the generic consequences of introducing B-C interaction leading 
to indirect positive feedback. Middle panels show the stripe (or lack of it) in the pure IFFMs and in the corresponding IFFM reinforced with the mutual- 
interaction. The dotted lines in the stripes represent attractor-switching, while the solid lines follow the attractor's movement. The Gray vertical bands 
qualitatively indicate the three regions, Low-High-Low, constituting the stripe in gene B. The yellow-background square panels are (C, B) phase plots 
generically corresponding to the three regions. There, the steady states (stable, S or unstable, U) are the intersections of nullcline curves (where one variable 
does not change in time - notation in the upper left corner). The black star indicates the initial condition close to the origin. The red arrows in these phase 
plots show that the nullclines move only horizontally or vertically in response to the morphogen gradient. When creating mutual inhibition, even a weak 
B-C inhibition allows single-stripe formation through attractor-switching process. When creating mutual activation, any value of B-C activation 
sharpens the stripe that already exists, and allows attractor-switching in the posterior border. The thin line passing through the unstable state shows the 
separatrbc, delimiting the basins of attraction for the two stable states. 



B-C inhibition are sufficient to allow for the multistability config- 
uration prone to stripe formation. 

In summary, the nullclines of the 2D simplified system act as 
dynamical building blocks for the generation of the stripe (see 
Supporting Information, Fig. S9 for a level-6 example). The strength 
of the interactions in the 3 -gene system determines the shape of the 
nullclines, ranging from a straight line to a S-shaped sigmoidal. The 
change in morphogen concentration translates into a simple shift of 
these nullclines, affecting the existence and accessibility of attractors. 
The topology of the network, through the sign of the interactions, 
determines the orientation of the nullclines and the direction of their 
shift. 

Discussion 

We have characterised analytically and numerically the interpreta- 
tion of morphogen gradient into a single-stripe of gene expression by 
network motifs based on two Incoherent Feed-Forward Motifs 
(IFFMs). Clear similarities were identified in the stripe-formation 



mechanism of the two IFFMs. In contrast, adding positive feedbacks 
to these two motifs generates different and counter-intuitive res- 
ponses. These differences can be summarised as follows. Firstly, 
the effect of the added interaction is independent of which positive 
feedback-loops already exist, i.e. the addition of these positive feed- 
backs is modular within a motif. By contrast, the same positive feed- 
back added to the two different motifs has different consequences on 
the stripe formation, i.e. the addition is non-modular between the 
motifs. 

Secondly, the addition of positive feedbacks is beneficial for the Ij 
motif (the beneficial-stripe motif), while it is neutral for the I 3 (the 
neutral-stripe motif). We traced the cause of neutrality to the 
dynamics conferred by the topology of the network motif. The 
phase-space reflects how the dynamics allows the Ij motif to main- 
tain the stripe for any strength of the added positive feedback. 
Conversely, the dynamics associated to the I 3 motif opens up new 
functional parameter regions, and thus provides parametric 
robustness. 
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Considering the IFFMs as single-stripe motifs, and the positive- 
feedbacks as border-sharpening motifs, the non-modular nature of 
their combination requires discussion. Recent studies have begun to 
hint towards the fragility of the dynamic-modularity assumption and 
the 'undesired' functional aspects of the cross-talks between motifs. 
For example, Ishihara & Shibata 38 show that combining (single- 
stripe) IFFMs into large networks gives rise to a new phenotype, 
the multi-stripe pattern. Complementary to the computational nat- 
ure of systems biology, the approach used by synthetic biology is 
empirical. Synthetic biology aims at building simple genetic circuits, 
and ultimately at assembling these well-characterised biobricks into 
intricate, customisable systems. Obtaining such reliable, complex 
gene circuits has proved to be more challenging than expected 39 . 
While modularity would be undeniably useful in synthetic biology, 
as it is in electronic circuits, it is still not clear to what extent complex 
biological networks can indeed be broken down into functional 
motifs. It has been argued that the existence or lack of such property 
might reflect different evolutionary conditions 40,41 . Nevertheless, the 
obstacles found in the straightforward scalability of modules reveal 
the necessity for further exploration of network dynamics, from 
simple motifs up to more complex networks. 

As the third and final result, the addition of the mutual-inhibition 
module stands out among all positive feedbacks by showing the 
largest expansion of the parameter space among all motifs studied. 
Moreover, these newly available parameter sets lead to high-express- 
ion single-stripe patterns with sharp borders, making them ideal for 
developmental patterning. Indeed, mutual-inhibition is a widely- 
employed pattern-formation motif in developmental gene net- 
works 610,12 . Other computational studies pin-point the numerous 
biological features the mutual-inhibition module brings about, 
especially the increase in sharpness and precision of pattern's bor- 
ders 22,26,4 " 5 . All these observations suggest that the mutual-inhibi- 
tion module truly defines a basic design principle. 

The IFFMs have been extensively studied for the variety of bio- 
logical functions they perform: time-dependent and dose-dependent 
biphasic response 7,46,47 , fold-change detection 48 , optimal information 
processing 49 and noise buffering 50 . Among the four possible IFFMs, 
architectures based on I] are the most abundant ones in bacteria and 
yeast 51 , a fact hypothesised to be due to its time-dependent charac- 
teristics of the IFFM. Additionally, far more I^based motifs than I 3 
have been reported for translating morphogen gradients into single- 
stripe expression patterns. As well-established examples we mention 
the triplets activin-Xbra-Gsc 15,26 , Bicoid-(knirps,Kriippel) -hunch- 
back 52 , Dorsal-(vein, rhomboid, sog)-Snail 53 , EGFR-Broad-Pointed 54 . 
We are drawn into speculating that this abundance might in fact be 
because the Ii motif is commonly forming part of a mutual-inhibi- 
tion motif, even if the second reressive interaction may have been 
missed in many cases. For the case of the activin-Xbra-Gsc, the 
mutual- inhibition has been indeed confirmed 26 . 

With respect to the general approach of our study, a few caveats 
are in order. For describing the integration of different transcription 
factors at promoter level we employed the sum & filter formalism. 
The choice of the formalism determine the types of topologies 
revealed by the subsequent analysis. For example, the I 2 motif 
(Fig. 2) could become a functional stripe-maker by changing to a 
constitutive promoter, an additional constant activator or alterna- 
tively a pronounced self- activation of the downstream genes. Under a 
different framework, such as filter & multiply 22 or a more general 
Hill-like formalism 24 , even the pure I 4 could show stripe formation 
due to the constitutive and AND-like features of the regulation func- 
tion. However, we chose to restrict to the current formalism given the 
wide-spread use of the sum & filter formalism in developmental 
biology, and in view of existing works on stripe-forming dynamical 
mechanisms 28 . 

Still within developmental biology, the biphasic dose-dependent 
behaviour addressed here under the term "single-stripe formation" 



relates to the positional-information mechanism. By this mech- 
anism, cells directly read and thus interpret morphogen concentra- 
tion-gradients leading to their subsequent differentiation in the early 
developmental stages. While there is now clear evidence that pattern- 
ing by morphogens is more complicated than previously thought 55,56 , 
our goal here was the understanding of minimal regulatory motifs 
with the aim of elucidating the mechanistic issues and underlying 
design principles of this patterning process. Experimental and com- 
putational studies alike continuously emphasise the extent to which 
the developmental program uses such functional motifs 6,21,57,58 . In 
other words, these observations support the hypothesis of design 
principles that conceives morphogenesis in terms of decomposable 
functional parts carefully assembled by evolution. 

However, a central question still remains, and that is why a 
specific network-design is chosen during evolution, when alterna- 
tive designs yield similar outcomes. The work of Cagatay et al. 59 is 
an iconic answer to this question, showing that the topology of a 
genetic regulatory circuit affects its noise characteristics, as net- 
works with identical deterministic outcomes can have profoundly 
different behaviours when noise is taken into account. In the cur- 
rent patterning context, Rodrigo & Elena 24 perform a first com- 
parison of stripe-forming IFFMs with respect to their noise 
response, showing that the four IFFMs are equally good at forming 
narrow stripes. Among all IFFMs reinforced with positive feed- 
backs, only the mutual-inhibition has been thoroughly studied 
and found to enhance the robustness of gene-expression bound- 
aries against variations in morphogen levels 44,45 . However, a 
detailed analysis and comparison of the noise response in IFFMs 
reinforced with positive feedbacks has not been performed yet. For 
discussion purposes, we performed a preliminary comparison 
between stripe-forming mutual-inhibition and mutual-activation 
modules. It reveals once more the robustness of the mutual-inhibi- 
tion module: while mutual-inhibition reliably maintains the stripe 
borders in spite of large fluctuations, mutual-activation shows a 
rapid shift or even loss of the posterior border (Supporting 
Information, Fig. S10). As suggested by other studies 60 , we believe 
this noise robustness may indeed be related to the parameter 
robustness revealed here in our study. In the light of all these 
observations, we consider a thorough stochastic analysis of all 
reinforced IFFMs as the logic and necessary future work that could 
explain the abundance of certain patterning designs in nature. 

Methods 

Simulation model. The parameters employed in the numerical simulations coincide 
with those from Cotterell & Sharpe 28 (electronic supplementary material): for the 
sigmoidal regulation function from Figure 2, we take a = 5, b = 1; the protein 
degradation rate 3 = 0.05; the morphogen gradient as the input into the A-gene 
follows an exponential increase as M — Id c , with I = 5,d = 0.982, and c, the integer cell 
number decreasing as c: 32 — 1. The promoters of the three genes are characterised by 
the same sigmoidal function. The steady states of the system are calculated from the 
ordinary differential equations: 

A= r — rrr\SA (!) 

l + exp(a-bW c ) v ' 



1 + exp[a - b(m AB A - m BB B - (o CB C)\ 



1 + exp[a — b{w AC A - m BC B — (OccC)\ 

The strength of the interactions are considered \coxy\ e [0> 10], and the initial 
condition, (A, B, C) = (0.1, 0.1, 0.1). The exhaustive search of stripe-generating 
parameters is performed under a sampling of A(o = 0.2. Sampling under a finer 
resolution maintains the results obtained, including the quantitative measures of the 
fraction of parametric volume (Supporting Information, Fig. Sll). We verified in 
such comparisons that there is no translation of the functional values, such that a ratio 
of 1 implies a complete superposition. Similarly, in a comparison between a level-4 
and a level-5 network, the common space to be compared is 4D. 

The definition of stripe in the numerical simulations considers the existence of a 
peak whose left and right border are set by a value of B Low = 10% of the maximum 
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value possible for the protein expression. This maximum expression is B Max — Vd, 
the inverse of the degradation rate. 
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